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Abstract 

As the concentration of intrinsic defects becomes sufficiently high in 0-deficient ZnO, interac- 
tions between defects lead to a significant reduction in their formation energies. We show that the 
formation of both O-vacancies and Zn-interstitials becomes significantly enhanced by a strong at- 
tractive interaction between them, making these defects an important source of n-type conductivity 
in ZnO. 

PACS numbers: 71.55.Gs, 72.80.Ey, 73.61.Ga 
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It is well-known that 0-deficient ZnO can easily become n-type even without the intro- 
duction of any intentional dopants. The mechanisms leading to the n-type behavior are, 
however, still controversial. Even though native defects resulting from thus deficiency have 
been excluded as the main source of the high free electron density, the n-type conductivity 
of ZnO is seen to be closely related to its 0-deficiency which manifests itself through the 
formation of 0-vacancy (Vq) and/or Zn-interstitial (Izn) defects. 

Theoretical first-principles studies of the formation enthalpies of point defects in ZnO 
have indicated that neither Vq nor Izn can lead to a high concentration of free carriers, 
since the most stable donor-like defect Vq has been shown to be a deep donor and the 
formation of the Izn shallow donor state has been known to be energetically far less favorable 
;han that of Vp when the Fermi level is close to the conduction band minimum (CBM) 
1 B. 0. a fl a a. 0„ the othe. ^a^a h d conta„.i„atio„ was proposed to be a. 

important cause of the natural n-doping [8|, 19|, [lOj • A potential problem with this proposal is 
that a high concentration of electron carriers is still observed even when H contamination is 
minimized or when H is removed , [isl. The sample annealed at 1100 °C in Zn vapor 



had an e^ 



cm ^ 



ectron concentration of 1.5x10^^ cm ^, while the untreated sample had 1.3x10^^ 



11| . Hydrogen easily diffuses out of ZnO at high temperature, and the H-related 



hyperfine structure observed in H-contaminated samples from the electron-nuclear-double- 

n n 

resonance study [Ij] was not observed in such heat-treated samples [15|. A meta-stable 
shallow donor state of Vq was suggested as an alternative source of the n-doping [l6|, but 
it is still controversial. 

As the concentration of defects becomes high, interactions between defects become un- 
avoidable and they need to be taken into account. It is known that donor-like and acceptor- 
like defects are attracted towards each other through their Coulombic interaction 1^]. How- 
ever, the defect-defect interactions between donor-like defects have been neglected in previ- 
ous theoretical studies of ZnO. In this Letter, we show that there can be a strong attractive 
interaction between two donor-like defects, the deep-donor Vq and the shallow-donor Izn- 
The driving force for the attractive interaction is the quantum mechanical hybridization 
between the electronic orbitals of their respective deep and shallow donor states, which at- 
tracts Izn toward Vq. The interaction significantly lowers the energy of the electronic donor 
orbital of Vq, as the distance between two defects decreases. The overall effect is a large 
reduction in the total energy of the system. As a result, the concentration of Izn can reach a 
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high enough level to explain the high concentration of electron carriers in 0-deficient ZnO, 
even if the Fermi level is close to the CBM. 

We performed density functional theory calculations, as implemented in the Vienna Ab 

examine interactions between 

defects, we employ a supercell geometry, in which two defects are located at various sites 
in the cell, and calculate the variation of the total-energy depending on their distance. We 
consider two low-energy donor-like defects: Vq and Izn- In n-type 0-deficient ZnO, the 
Fermi level is close to the CBM, and the neutral charge state of Vq is stable, since Vq is 
a deep-donor defect, but the shallow-donor Izn is stable in a (2-|-)-charge state. Therefore, 
only interactions between these charge states are focused on in this work. The results of our 
calculations are shown in Fig. [H where the variations of formation enthalpies as 
a function of defect-defect [a-a') separation (r) are shown. The Fermi level is assumed to 
be located at the Vq deep donor level. 

We first consider interactions between two Vq's and, separately, between two Izn's. It is 
found that the total energy increases as the distance between the defects decreases, indicating 
a repulsive interaction between same-type defects. As shown in Fig. [1], the energy of two 
Izn's increases sharply as they get closer to each other, whereas that of two Vq's increases 
slightly. The repulsive interaction between Izn's is due to the expected Coulombic repulsion. 
The variation of the energy is fit well using the screened Coulombic interaction functional, 
■ ^ — , where the fitting parameters are the charge state Z=2, the dielectric constant 
e=8.65, and the screening length A=1.13 A. The weak repulsive interaction between Vq's is 
related to a strain effect. The inward- relaxation of neutral Vq induces a tensile-strain into 
the surrounding lattice. In the paired state, the strain does not become fully relaxed. 

An interesting finding is that there is a strong attractive interaction between Vq and Izn- 
It is surprising that an attractive interaction is induced between two donor-like defects. The 
formation enthalpy of Vo-Izn pairs as a function of their separation is shown in Fig. [H The 
total-energy decreases significantly, as the separation becomes smaller, indicating that Izn 
is stable with respect to a coupling to an on-site Zn atom surrounding Vq. In Fig. [2]^a) 
and (b), the possible sites for Izn around Vq are indicated; the octahedral interstitial sites 
(01-08), the tetrahedral interstitial sites (T1-T3), and the Vq site (V). The most stable site 
is the 03 octahedral interstitial site for which the atomic structure is shown in Fig. M^c). 
The binding energy is 0.52 eV. There are three equivalent 03 sites around a Vq. When Izn 
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is trapped just at the oxygen vacancy site V (as an antisite defect), the 0.44 eV binding 
energy is shghtly smaller. 

In order to understand the driving force of the attraction between Vq and Izn, we exam- 
ined the change of the electronic structure, depending on distance. The electronic structure 
of the coupled state is significantly modified as the separation between the two defects be- 
comes smaller. It is noted that the defect levels of the Vq (both the occupied ai gap state 
and the unoccupied t2 state) are lowered while those of Izn state are raised. The calculated 
results for the electronic structure are described by the local density-of-states from the four 
Zn atoms around Vq and from the Izn atom, as shown in Fig. [3](a)-(e). Figure [3](g) de- 
scribes the lowering of the Vq state as a function of the distance. The lowering is noted to 
be quantitatively similar to the change of the formation enthalpy in Fig. [TJ This indicates 
that there is a strong interaction between these two defect-levels. 

The driving force for the attraction originates from the hybridization between the deep 
donor state and the shallow donor state. The change in electronic structure can be simply 
depicted as in Fig. [3t^f). The defect level of the shallow donor Izn is located inside the 
conduction bands and the Vq level is located deep within the band gap. Both defect orbitals 
mainly come from the Zn-4s orbitals. As the distance between the two defects becomes 
smaller, a hybridization between the two defect orbitals is induced, as a result of which the 
lower-lying Vq level is lowered, and the upper Izn level is pushed-up. Since the Vq level 
is occupied by two electrons while the Izn level is empty, the hybridization lowers the total 
energy. As the defect-defect separation becomes smaller, the hybridization becomes stronger. 
This gives rise to an attractive interaction between the deep and shallow donor states. As 
shown in Fig. [3](g), the change of the Vq level is well- fit by an exponential function: doc'^'^""* , 
with the fitting parameters, (5o=0.47 eV, which is 5 at r=0 A, and a*=4.29 A. The wave 
functions of the defect levels of the donor-like Vq and Izn can be approximately written by 
g-r/av g^j^^ e"''/"', respectively, where ay and ai are the effective Bohr radii of the defect 
electron densities. The deep level of Vq is localized nearly within the second oxygen nearest 
neighbor, av~5.5 A, and the shallow donor state is delocalized as ai=aoe/m*=20 A within 
the simple effective mass theory, where the Bohr radius ao=0.53 A, the dielectric constant 
e=8.65, and the effective mass of electron m*=0.23. The fitting value of a* can be obtained 
from the product of the two wave functions, {ijjy\i(ji), with av=5.5 A and ai=20 A, i.e., 
l/a*=l/av+l/ai=l/(4.3 A). There is just a slight deviation between the lowering of the 
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total energy (Fig. [T]) and the Vq level [Fig. M^g)] from the fitting curves, which is related 
to the local strain effect. The closest Vo-Izn pah (V) shows the largest level drop 0.403 eV 
[see Fig. El^e) and (g)], and the lowest energy configuration is found at the 03 site, because 
the 03 interstitial volume is more spacious as it is located behind the inward-relaxed Vq 
[see Fig. [2](c)], and thus the accommodation of Izn is easier. The interaction via orbital 
hybridization between the defects is more significant than the elastic interaction. 

The attractive interaction between Vq and Izn can lead to the coexistence state of both 
Vq and Izn rather than the presence of only Vq in 0-deficient ZnO. We estimated the con- 
centrations of Vo-Izn pair, relatively to isolated defects, and also the Fermi level through 
the calculations of the formation enthalpies of defects under the thermal equilibrium condi- 
tion. The average formation enthalpy of a defect a under the interaction with a'-defects is 
calculated by 

^a=^°+5Z$^^a,«'W, (1) 

r a' 

where is the formation enthalpy in the isolated state without interaction, and Ua,a' {r) is 
the interaction energy between a- and a'-defects separated by r. In Fig. [H the calculated 
= + + Ua,a'{r) are shown, by which Ua,a'{r) can be estimated. For the 
numerical simplicity, particularly in the limit of extreme high concentrations of defects, 
we adopt a mean-field concept, and Qa can be approximated by + J2a' ^a,a', where 
Ua,a' = / dVUa,a'{f)na' is the mean interaction energy between defect a and a', and f/a,a'(r) 
is approximated by the fitting curves drawn in Fig. [T] Here ?7,„/ is the density of a defect 
a'. When we employed the calculated values for Ua^a'i'r') shown in Fig. [H the overall results 



show little change. Now we follow conventional method (20| to determine the concentration 
of each defects by Ua = Noe~^^°' , the Fermi level (^f), and the carrier concentration [rie) 
by charge neutrality condition. Nq is the number of available sites for defect formation: 
about 4.22x10^^ cm~^ in ZnO. The concentration of electron carrier depending on ep is 
determined by the effective density of electronic states of the CBM (Nc) of 2.8x10^^ cm~^ 
at room temperature and a conduction-band electron effective mass of m*=0.23. Since the 
concentration of hole carriers is minor, so the results are insensitive to the density of the 



2lj. The 



state around the valance band maximum, we used an effective mass of m^=0.5 
formation enthalpies of the isolated Vq and are 1.68 eV and 2£f-0.91 eV, respectively, 
under the Zn-rich condition, where the Zn chemical potential is chosen to be the energy of 
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a Zn atom in the Zn hep metal 19|. The crystal growth condition varied between 0-rich 
and Zn-rich conditions is represented by the O chemical potential, which affects the defect 
concentration. 

In Fig. m the concentration of Izn (ni) and the Fermi level are shown with respect to the 
concentration of Vq (ny) at a growth temperature of 900 K. When we neglect the interactions 
between defects, the Fermi level is calculated to be pinned at a mid-level between the CBM 
and the Vq level, and ni is limited to much less than 10^^ cm~^, as conventionally indicated. 
When we include the interactions between defects, the result is significantly changed. Under 
0-deficient condition, ni can become close to riy, i.e., the coexistence of the Vq and Izn in 
the 0-deficient ZnO becomes remarkable. Since Izn is a shallow-donor defect, Fermi level can 
be raised to levels much higher than those expected in the presence of only Vq's and thus 
the 0-deficient ZnO can be heavily n-type doped. These findings show how the 0-deficiency 
itself can be the source of the n-type conductivity in ZnO. When ny rises up approximately 
10^° cm~^, where the average distance between Vq's is about 27 A, ni can reach about 10^^ 
cm~^, and the Fermi level approaches the CBM with an electron carrier concentration rig of 
2x10^^ cm-3. 

Recently, Vlasenko and Watkins have reported the existence of a complex defect consisting 
of a deep donor Vq and a shallow donor in ZnO through the optical detection of electron 
paramagnetic resonance (EPR) in electron- irradiated ZnO 2^. They identified an EPR line 
whose g value is accurately given by the average of the g values of the EPR signals from Vq 
and a shallow donor. The coincidence was suggested to be the result expected for a closely 
located defect pair, i.e., Vq and a shallow-donor. The pair may be explained by the Vo-Izn 
pair reported in this present study. It should be also noted that the defect level of Vq can 
be deeper below the midgap of ZnO when coupled with Izn than isolated. 

In conclusion, we find a strong attractive interaction between Vq and Izn, whose driving 
force originates mainly from the hybridization between the orbitals of the Vq deep and the 
Izn shallow donor states. An Izn defect preferentially forms around a Vq defect as a result 
of the attractive interaction. We find that in seriously 0-deficient ZnO, the concentrations 
of Vq and Izn can be comparable. Therefore, 0-deficiency itself can be the source of the 
ri- type conductivity in ZnO. 
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FIG. 1: (color online). Calculated formation enthalpies of the defect pairs as a function of the 
distance. Attractive interaction of the Vo-Izn pair (^)) and repulsive interactions between the 
same types of defects, Vq-Vq (H) and Izn-Izn (^) are shown. 
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FIG. 2: (color online). Available atomic sites for an Izn near a Vq are indicated in the ZnO 
wurtzite structure shown from (a) the top (along the polar axis) and (b) from the side view, (c) 
Atomic structure of the most stable Vq-Izh pah (03). The large and small atoms indicate the Zn 
and O atoms, respectively. 
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FIG. 3: (color online). Calculated local density-of-states of electrons near Vq (red) and near Izn 
(blue) for (a) the separated Vq and Izn (a), and for the Vq-Izh pairs of (b) 07, (c) 05, (d) 03, 
and (e) V of which sites are indicated in Fig. [2][a) and (b). (f) Schematic diagram of the defect 
levels of Vo, Izn, and Vo-Izn- Vq and I^^ denote the Vq and Izn in Vo-Izn pair, respectively, (g) 
Calculated level drops 6 of the Vq state as a function of the distance. 
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FIG. 4: (color online), (a) Fermi level and (b) equilibrium concentrations of Izn calculated by 
considering the defect interactions (red) and neglecting the interactions (blue) as a function of Vq 
concentrations. Theoretical values of the CBM and the Vq level are indicated in (a). Equilibrium 
electron carrier concentrations {n^ are also plotted in (b). 



12 



